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Abstract 

This paper applies the analytic forms of a recent non-perturbative, manifestly gauge- and 
Lorentz-invariant description (of the exchange of all possible virtual gluons between quarks (Q) 
and/or anti-quarks (Q) in a quenched, eikonal approximation) to extract analytic forms for the 
binding potentials generating a model Q-Q "pion" , and a model QQQ "nucleon". Other, more 
complicated Q, Q contributions to such color-singlet states may also be identified analytically. 

An elementary minimization technique, relevant to the ground states of such bound sys- 
tems, is adopted to approximate the solutions to a more proper, but far more complicated 
Schroedinger/Dirac equation; the existence of possible contributions to the pion and nucleon masses 
due to spin, angular momentum, and "deformation" degrees of freedom is noted but not pursued. 
Neglecting electromagnetic and weak interactions, this analysis illustrates how the one new pa- 
rameter making its appearance in this exact, realistic formalism may be evaluated, along with a 
qualitative estimate of the lowest quark mass and the energy contributions to the ground states of 
the gluon fields, from a knowledge of the pion and (ground state) nucleon masses. 
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I. INTRODUCTION 

This paper should be viewed as a continuation of that of Ref. in which a new formalism 
for analytic QCD calculations was suggested, and which can serve to identify an effective 
potential as the result of the exchange of multiple, virtual gluons between Q's and/or Q's. 
The word 'effective' here represents the type of potential one would expect to find were the 
scattering pair of conventional Abelian particles, moving at large relative velocities in their 
CM, as appropriate in an eikonal context. And the word 'Abelian' is most relevant in this 
QCD situation, for the analysis of Ref. l| clearly shows that for large impact parameters 
the scattering of such non- Abelian Qs and/or Qs is effectively coherent; only as the impact 
parameter decreases do progressively stronger color fluctuations appear and decrease the 
the effective coupling, which then becomes the statement of non-perturbative asymptotic 
freedom in the variable conjugate to impact parameter. 

In the present paper, we shall be interested in binding, or "restoring" potentials, effective 
potentials which prevent bound particles from remaining at relatively large separations. 
Since large impact parameters are relevant to large, three-dimensional distances, what is 
needed is, firstly, a statement of the relevant formulae of Ref. [l] at large impact parameter; 
and, secondly, the corresponding eikonal function found in that limit. It is then a relatively 
simple matter to determine an effective potential which would produce the same eikonal. 

There is one truly profound statement which deserves careful mention, as in Ref. [l| 
and in the paper immediately preceding this one [3], which is that the quanta of the the 
non-Abelian QCD fields differ from quanta of Abelian fields in that the latter quanta, in 
principle, satisfy standard quantum-mechanical measurement possibilities: position can be 
specified exactly at the expense of momenta, and conversely. Quarks and gluons, in contrast, 
are bound objects; and in the context of Q/Q scattering, while their longitudinal and time 
coordinates (or their longitudinal momentum and energy) can be estimated in terms of 
that of the hadron into which they are bound, their transverse momenta, or the (impact 
parameter) distance between them, can never be measured with precision. It is precisely 
here that one is forced to bridge the gap from "ideal" to "realistic" QCD; and until a more 
fundamental field theory of QCD is developed, which contains such imprecision from its very 
beginning, one will be forced to make a "phenomenological" transition into the "realistic" 



world of physical QCD, as in Ref. {3] 



The need for such a transition was made quite clear in Ref . |l| , where one of a variety of 
forms was suggested as appropriate, when it became clear that a 5^ (b) must be replaced by 
a smoothly-falling function of impact parameter, e.g., ip(b) ~ e - ^ 2 , although no definite 
choice was needed to reach the conclusions of that paper. In the present context of binding, 
a specific choice is made; but it should be emphasized that this choice is phenomenological, 
if reasonable: it works, generating a binding potential V(r) proportional to £ -r 1+ ^, where £ 
is a small, real and positive parameter. The appearance of £ was unexpected, and represents 
a new quantity which should appear naturally in a more profound version of QCD; it is here 
specified by a minimization procedure of the ground-state's energy level of our model pion. 

It should also be emphasized that the transition to "realistic" QCD has nothing to do 
with the quenched, eikonal approximations used to give the simplest possible, non-trivial 
example of our new approach to gauge-invariant, non-perturbative QCD. Exactly the same 
necessity of such a transition may be seen in the second paper of this series j^j], where the 
forms of the results displayed in Ref. [ij are showed to hold without approximation and 
without exception. 

Our present discussion will follow that of Ref. [l| in extracting, under quenched ap- 
proximation, the eikonal and a corresponding potential function for the equal-mass Q-Q 
contribution to the model singlet state we here call the "pion". To make a similar statement 
for the QQQ contribution to the model singlet ground state we call the "nucleon", it will 



be useful to recall the eikonal forms appropriate to three-body scattering amplitudes 



3,3, 



while for application to higher Q/Q contributions to these and other hadrons, the discussion 
of N-body eikonals of Ref. % Iq] will suffice. 

We refer the reader to the unique method employed in Refs. [l, 2] to convert Schwinger's 
conventional, gauge-dependent functional representation for the QCD generating functional 
into an exact, related representation, in which all processes corresponding to Q and Q 
interactions by the (non-perturbative) exchange of all possible gluons - including Q and Q, 
cubic and quartic interactions - are rigorously and manifestly gauge invariant. Attention is 
then drawn to the quenched Q-Q scattering amplitude in eikonal approximation, and to the 
manipulations and (strong-coupling) simplification used to obtain the integral of Eq. (44) of 
Ref. JJ, and the replacement of its integrand by Eq. (60). We begin our extraction of the 
Q-Q binding potential from this point. 



II. QUARK BINDING POTENTIAL 



In this Section, we shall be interested in binding, or "restoring" potentials, effective 
potentials which prevent bound particles from remaining at relatively large separations. 
Since large impact parameters are relevant to large, three-dimensional distances, what is 
needed is, firstly, a statement of the relevant formulae of Ref. l| at large impact parameter; 
and, secondly, the corresponding eikonal function found in that limit. It is then a relatively 
simple matter to determine an effective potential which would produce the same eikonal. 



A two-body eikonal scattering amplitude will have the form 

is 
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where s and t denote the standard Mandelstam variables, s = — (P1+P2), t = —(jPi—Pi) 
q 2 in the CM of QQ, and X(s, b) is the eikonal function. Following the arguments of Ref. 
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(2) 



one reaches 

e iX = M f d"x30 [detGtf • x)" 1 ] 1 • & ■ exp 
and the normalization is such that for ip(b) — > 0, exp [zX] — > 1. All the integrals over 

d n ai [ d n a R [ d n fii / d n n n (3) 



are properly normalized, and they connect the Qf, dependence of Eq. §Z§ with the Af , 
Aj Gell-Mann matrices placed between initial and final state vectors, which define the type 
of amplitude desired, Q-Q scattering or Q-Q scattering, etc. 

As in Ref. we arrange the 8-dimensional integral over J d n Xo3 into the form of a radial 
integration over the magnitude of X03, now called R, and normalized angular integrations. 
We suppress the latter, along with all normalized j d n Qi j d n Qi dependence, and rewrite 
Eq. (j2J) in the form 

e«(*> 6 )=W dRR 6 - 5 e l ^ +l ^ (4) 
Jo 

1 00 • R 2 

where the normalization of the i?-integral is given by AT'' 1 = J™dRR 6 - 5 e i£! t, and where 
(g) is a shorthand denoting the effective coupling AFTER all angular and color integrations 
and matrix elements have been calculated. In this calculation, only the linear dependence 
of gip is needed, while the non-vanishing matrix elements between singlet color states will 
define Q-Q scattering. 



As noted in Section [J we wish to compare this eikonal function at fairly large impact 
parameter with that corresponding to the scattering of two, interacting Abelian particles, 
and then infer what form of effective potential corresponds to that eikonal. We therefore 
assume the usual formula [8] relating an eikonal to a specified potential, 

/oo 
dzV(b + P L z) = 7 (s)X(6), (5) 
-oo 

where j(s) depends upon the CM energy of the scattering particles (and the nature of their 
interaction), and Pl is a unit vector in the direction of longitudinal motion. Eq. (jHJ) is true 
for non-relativistic and relativistic scattering, in potential theory and in field theory; for 
large b, our eikonal will turn out to be only weakly dependent on s, and the ^-independent 
contributions do not contribute to V(r). To obtain V(r) from a given X(6) is then simply 
a matter of reversing the usual calculation: one computes the Fourier transform X(&j_) of 
X(6), and then extends k± to \k\ = [k 2 ± + kj]*. The three-dimensional Fourier transform of 
X(fc) is then proportional to V(r). Following the discussion of Ref. jsj], we assume that X(6) 
- which depends on the assumed form of ip(b) - is a function of b 2 ; and although we begin 
by asking for the form of ip(b) for large b, it can be shown that the small-6 contributions to 
the Fourier transforms are unimportant for the large-r behavior of V(r). 

There is one further point which requires consideration, for the above remarks are valid 
when the process involves only the scattering of two particles. But when initial energies are 
high enough such that inelastic particle production occurs, then the potential corresponding 
to such production must have a negative imaginary component (so that the corresponding 
S-matrix, exp [— iHt], will have a decaying time dependence), thereby conserving probability 
as a diminution of the final, two-particle state when production becomes possible. And this 
is also true when the two initial particles have the possibility of binding, and forming a state 
which was not initially present; the probability of the final scattering state must diminish. 
In other words, the eikonal function one expects to find may always be characterized by a 
complex potential of form Vr — iVi, where in the present case, V\ is that potential which can 
bind the Q and Q into a pion; and that potential, in a scattering calculation, must appear 
with a multiplicative factor of —%. Symbolically, zX — > i{V^ — iV\), and for this binding 
calculation, zX — >• Vb. 

Return to Eq. ([1]) and assume q ^ 0, so that the " 1" term of the integrand contributes a 
5^ (q) which vanishes. This is a physical assumption: Q and Q are not bound rigidly, with a 



unchanging impact parameter; rather, there is a constant, if small, "in-and-out" transverse 
motion, which represents the bound state; and hence q ^ 0. Then, 

T ~ J d 2 b e^V x(b) (6) 

where exp [iX(fe)] is given by the normalized Eq. (j3J). In Ref. [l| we imagined that (fib) would 
be a smooth function which vanishes for large b, such as ~ exp [— (fib) 2 ], but it turns out, 
below, that an appropriate choice is 

(fib) = <p{0) e~^ 2+ \ (7) 



where £ is real, positive, and small. In the preceding paper, Ref. [3(, a detailed argument 
for such a "realistic" choice of ip(b) was given, leading to the generic Eq. (21) of that paper. 
For any such choice of positive £, ip becomes small when fib ^> 1. 

For large b it is then sensible to expand the exp [i(g)(p(b)/R] term of Eq. (J4]), retaining 
only the linear (g)<p dependence. This gives, in place of Eq. (j3J, 



e «W ~ j\f' / di?i? 6 V^ 
= l + iKAf'(g)ip{b) + - 



(8) 



where 



kN' = N' dRR^e^ (9) 
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and M' = 2(4i) 4 /r(^). Remembering that both sides of Eq. (jHJ) are to be integrated over 
J d 2 be l ^' b , for q ^ we can associate 



3 iX(6) 



iKAT{g)ip{b), (10) 



or 

iX(6)=ln[^(6)] + ..., (11) 

where In [yj] is large and y> is (effectively) small. Choosing (p(b) = ip(0)e~^ 2+ , we then 
obtain 

iX(b) = -(fib) 2+ t + ■ ■ ■ (12) 



It is then convenient to use the integrals [10] 

f°° 1 // 1 

J dxx^J (ax) = 2^a- 1 -^T(- + ^)/T(- 



2' 



H < 



2' 



and 



dxx M 1 sin(x) = T(fi) sin(^), |Re(/x)| < 1, 



(13) 



(14) 



in which, except for obvious poles, the Gamma functions are analytic in /i, and can be 
continued to the needed values. With 



iX{k ± ) = - / d 2 &(/^) 2+ «e iA ^ 



(15) 



/■oo 

-(2tt) / dbJ (k ± b)/i 2+ tb 3+ t, 
Jo 



and the use of the doubling formula for Gamma functions, working everything through, one 
finds 



V B (r) = -'— r 1 * T ^_i } r(-2 - sin(^), (16) 



and for small £, £ <C 1, 



V B (r)~^(iir) 1+ Z, (17) 

which may be compared to the results of several machine groups. There remains the question 
of what values should be assigned to £, and to the unknown parameter /i/wiQ; and these 
will be dealt with in the next Section. 

To calculate the corresponding, effective, restoring potential when one of the three quarks 
contributing to the QQQ bound state is separated from the other two, or, more simply, when 
all three quarks are forced to separate from each other, one may refer to Eq. (22) of Ref. Q], 
giving the eikonal amplitude corresponding to Coulomb three-particle scattering. Here, the 
specifically Coulomb parts of this amplitude may be replaced by the single Halpern integral 
which connects the three quarks to each other, as it connected the Q and Q to each other in 
the calculation above. If particle 2 of this formula enters the scattering with zero transverse 
momentum in the rest frame of particles 1 and 3, the amplitude simplifies to 



7? ik ~ / d% 2 I d% 2 e i ^ 632+i ^ h2 



<5>( y b 12 ,b 32 ,b 13 ) + # s 



(18) 



where 613 = 612 — 632, irrelevant kinematic factors multiplying these integrals have been 
suppressed, and ^3 denotes the combination 



1 + 



$(612, 00, 00) 



+ 



1 - $(00,632,00) 



+ 



1 — $(00, 00, b 



13 



(19) 
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so that 7^ eik is a completely "connected" amplitude. 

Writing $ as exp [zX(6i2, 632)], its defining integral may be rewritten as 

poo 

Af' / dRR 6 5 exp 
Jo 

and one may ask for the form Eq. fl20|) takes when one or more of the b^ becomes large. 

When any one of the 6^ becomes large, its <f(bij) becomes small, and its exponential term 
may be expanded. The simplest situation is when they all become large, so that an expansion 
of Eq. (|20l) is relevant. The first non-zero and leading term in such an expansion which 
contains the (p dependence, and hence the ^-dependence, of all three quarks contributing 
to the singlet QQQ state is that term in which the expansion yields 

1 + kM' \{g) l2 ip{b l2 ) + (g) 3 Mh 2 ) + (g) l3 <p(h 3 )] + • • • , (21) 

so that, suppressing all normalized angular and color integrations, for q% ^ 0, q 3 ^ 0, 

e iX(fel2 ' fe32) * kAT MnVibn) + (9)^32) + (g) 1 M b ^)} + • • • ■ (22) 

If particles 1 and 3 are now additionally separated, keeping the distances 612 and 632 essen- 
tially fixed, then 

e iX (bl2) M ^ Kfir^tpQto) + ■■■ , (23) 
iX(b 12 ,b 32 ) ^ In [<p(b 13 )] H , 

and we are effectively back in the "pion" situation, where the large- impact parameter eikonal 
of Q-Q scattering was calculated. Clearly, Eq. fl23|) suggests that when any two of the three 
quarks which contribute to the singlet "nucleon" are well separated, there results a restoring 
potential V(ry-) of the same form as that of Eq. ( TlTl) . There are, of course, corrections to 
this confining potential, those which are obvious from the approximations used above, as 
well as those corresponding to spin and angular momentum, which have been completely 
neglected. 

III. ESTIMATION OF THE "MODEL PION" MASS 

This Section contains an estimate of the effects of the above binding potential for the 
simplest case of the "pion", a ground state of the Q-Q system calculated using the simplest 
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R 

+ i [(g)iMbi2) + {g) 3 Mh 2 ) + (g) 13 ip(b 13 )] 



(20) 



minimization technique as well as a second minimization with respect to a ratio of terms 
introduced in this analysis. In this simplest, non-relativistic estimation, the Hamiltonian of 
this system is given by 

H = 2m+— p 2 + V(r), (24) 
m 

where 1/m denotes the inverse of the reduced mass of this equal-mass quark system. The 
ground-state energy is then given by the replacement of p by 1/r, and the subsequent 
minimization of the eigenvalue of this Hamiltonian with respect to r. Assuming £ <C 1, one 
obtains 

-(IS 1 ". 

and 

1 /3 

E = 2m + ifr (- !L) , (26) 
2 \£ in/ 

Here, \x and £ are the two parameters required by this analysis in order to have a non-zero 
binding potential; and while their existence has been introduced as a necessary assumption 
for the interacting Q-Q system, there is another unknown quantity, the quark mass m, which 
must play an important role. If we take the position that /i and £ will appear as a result 
of a more fundamental QFT (in which the transverse momenta or position coordinates of 
the quanta of the QCD fields can never be measured with precision), then we are free to 
consider which values of m, or of the ratio fi/m, can lead to the lowest value of E. In so 
doing, we are asking if there might exist a dynamical reason for the relatively small value of 
the pion mass, in this approximation to the actual pion; that is, whether approximate chiral 
symmetry has a dynamical basis. 

Following this approach, and using the small-£ limit of Eq. (1261) , 

-m/ V 2 / \m/ 



E = /I 



(27) 



upon variation with respect to x — [i/m, E (x) will have an extremum at 
and since 



d 2 E n (x) 



dx 2 



■I'D 



( |) 3/2 (2-|) >0, (28) 

-1/2 



that extremum is a minimum. Substituting the value of [ifm = 2 3 / 4 f |J into Eq. (127)1 
yields 

J E = / i£ 1 /2 2 -i/4) [1 + 3]) (29) 



from which one infers that there is three times as much energy in the gluon field and Q- 
Q-kinetic energies as in the quark rest masses. Intuitively, one expects that E ~ m n ~ /i, 
which suggests that £ ~ v2/16; and, finally, one then has an estimate of the quark mass, 
in terms of £ and fi. The value of m n may be inferred from the fall-off of the (approxi- 
mate) Yukawa nucleon-nucleon potential for moderate impact parameters, corresponding to 
separation distances in the " one-pion-exchange" region of that potential; and this will be 
discussed in a subsequent paper. 

Of course, this double-minimization calculation of Eq cannot be exactly identified with 
the precise pion mass, because of the approximations made above, as well as the omission of 
more complicated singlet terms, such as the contributions coming from QQQQ terms. But 
the Physics seems to be reasonably correct. 

A similar minimization analysis can be made for the nucleon ground state, but there then 
are other degrees of freedom which should be included, and treated properly by a serious 
attempt at a solution of such a three-body, relativistic problem. This we must leave to 
others, whose numerical abilities far exceed our own. 

IV. SUMMARY 

In this paper we have indicated how analytic binding potentials "naturally" appear in 
terms of the eikonal amplitudes considered, as the simplest, non-trivial examples of the new, 
gauge-invariant approach to analytic QCD. In the necessary shift from "ideal" to "realistic" 
QCD, we have made use of a simple formalism [3] which incorporates transverse imprecision 
in terms of a small, real, and positive parameter £, whose numerical value, and that of the 
ratio of /i/m can be inferred by a comparison with the measured pion mass. Even though 
£ appears small, the contribution of the binding potential to the lowest bound-state-model- 
pion mass is large. One can expect the same sort of relationship between the nucleon mass 
and the parameters £, fi, and m to occur for the QQQ model of a nucleon. 

The next paper of this Series will deal with Hadronic Scattering Amplitudes, which em- 
ploys the same set of parameters, and which will produce effective, Yukawa-type exponentially- 
decreasing scattering potentials between bound, singlet combinations of quarks and anti- 
quarks, scattering at distances larger than typical impact-parameter distances of the binding 
calculation of this paper; and the exponential fall-off of the scattering potentials will also 
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depend on the same set of parameters. To this end it will be necessary to enlarge the 
scope of the non-perturbative terms included by removing, in a simple, illustrative manner, 
the quenched approximation of the binding calculation above. We emphasize that these 
calculations incorporate the exchange of all possible gluons between realistic quarks and 
anti-quarks in sums over all possible, relevant Feynman graphs: the functional techniques 
used here are intuitive, simple, and manifestly gauge- invariant. 
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